Quantum model of microcavity intersubband electroluminescent devices 
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We present a quantum theoretical analysis of the electroluminescence from an intersubband tran- 
sition of a quantum well structure embedded in a planar microcavity. By using a cluster factorization 
method, we have derived a closed set of dynamical equations for the quantum well carrier and cavity 
photon occupation numbers, the correlation between the cavity field and the intersubband polariza- 
tion, as well as polarization-polarization contributions. In order to model the electrical excitation, 
we have considered electron population tunneling from an injector and into an extractor contact. 
The tunneling rates have been obtained by considering the bare electronic states in the quantum 
well and the limit of validity of this approximation (broad-band injection) are discussed in detail. 
We apply the present quantum model to provide a comprehensive description of the electronic 
transport and optical properties of an intersubband microcavity light emitting diode, accounting 
for non-radiative carrier relaxation and Pauli blocking. We study the enhancement of the electrolu- 
minescence quantum efficiency passing from the weak to the strong polariton coupling regime and 
compare it with the free-space case. 



In the last two decades, the fundamental research 
on the physics of intersubband transitions in semicon- 
ductor quantum wells has enjoyed a considerable suc- 
cess and also led to novel applications in quantum 
optoelectronics [l|. Recently, reflectivity experiments 
[1, S S] have demonstrated that by embedding a doped 
quantum well structure in a planar microcavity, it is 
possible to achieve the strong coupling regime between 
an intersubband transition and a cavity photon mode, 
provided that a dense enough two-dimensional electron 
gas populates the fundamental quantum well subband. 
The interaction between a bright intersubband excita- 
tion and a cavity photon is quantified by the so-called 
vacuum Rabi frequency. The strong coupling regime oc- 
curs when the vacuum Rabi frequency exceeds the elec- 
tronic and photonic losses. In such a regime, the normal 
modes of the system are cavity polaritons, half-photon 
half-intersubband excitations. In this kind of system, it 
is even possible to reach an unconventional ultra-strong 
coupling regime, i.e. a vacuum Rabi frequency compara- 
ble to the intersubband transition frequency[5|, |6|, |7| . 

The interplay between judiciously quantum engineered 
intersubband transitions and vertical electron transport 
is the essence of the so-called quantum cascade electro- 
luminescent devices and lasers, which are unipolar opto- 
electronic sources emitting in the mid and far infrared 
portion of the electromagnetic spectrum [^, 0, Ufll- A 
new kind of microcavity-embedded quantum cascade de- 
vices in the strong coupling regime was proposed in Ref. 
I 111 . The first experimental demonstrations of a microcav- 
ity quantum cascade photovoltaic detector [12] and of an 
electroluminescent device in the strong coupling regime 
have been recently reportedp^. 
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This promising research topic is in its very infancy 
and many interesting theoretical questions need to be 
addressed. In Ref. [g, intersubband polariton electrolu- 
minescence has been analytically treated within a simpli- 
fied Hamiltonian model based on the following assump- 
tions: (i) only the bright intersubband excitations have 
been taken into account, while dark excitations have been 
neglected; (ii) only the low excitation regime has been 
considered, in which the bright intersubband excitations 
have been approximated as bosons; (iii) the electronic 
coupling to the intersubband polarization field has been 
modeled through a phenomenological reservoir of bosonic 
excitations. In this work, we will attempt to treat the 
same problem starting directly from the fermionic Hamil- 
tonian for the quantum well carriers. This approach can 
give us useful insight to understand the physics obtained 
by relaxing the assumptions used in Ref. |]g and to grasp 
which intrinsic factors ultimately determine the quantum 
efficiency of these strong coupling emitters. On one hand, 
the large values of the vacuum Rabi frequency could in- 
duce a very fast and efficient emission of photons. On the 
other hand, the large density of dark intersubband exci- 
tations created by the injection current and the Pauli 
blocking in the densely populated fundamental subband 
could suppress such enhancement. 

We would like to point out that from a theoretical 
point of view, a description of the considered system in 
terms of the fermionic carrier operators makes the system 
Hilbert space much larger than within a bosonic model. 
In this paper, we have followed an approach based on a 
truncation of the infinite hierarchy of dynamical equa- 
tions for the operator expectation values, allowing us 
to describe many relevant aspects of the intersubband 
microcavity electroluminescence. However, some of the 
non-perturbative features obtained analytically within a 
bosonic model [6| can not be accounted for within the 
present treatment. Different fermionic approaches based 



on exact diagonalization methods [lj| are eventually nec- 
essary for further refinements. 

In this paper, we present a quantum model of the spon- 
taneous photon emission from an electrically-excited in- 
tersubband transition of a quantum well structure em- 
bedded in a planar microcavity mode. Here, we will 
consider the case of an incoherent electron transport, 
where the quantum well electron populations in the two 
subbands have a tunneling coupling to an electronic in- 
jector and to an extractor. The tunneling rates have 
been obtained by considering the bare electron states in- 
side the quantum well. The domain of validity of this 
approximation will be discussed in detail. The present 
theoretical model is applied to describe the incoherent 
electron transport and electroluminescence of an inter- 
subband microcavity light emitting diode in the strong 
coupling regime. The paper is structured as follows. In 
Sec. m we describe the system and introduce the second 
quantization Hamiltonian describing electrons in the two 
conduction subbands and photons in the fundamental mi- 
crocavity mode. In Sec. [Ill we present a closed set of dy- 
namical equations for the one-time expectation values of 
operator products, describing photon and carrier popu- 
lations as well as intersubband polarization-polarization 
and polarization-field correlations. These equations have 
been obtained through a cluster expansion, whose de- 
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tails are reported in Appendix [XI In Sec. IIIIl we discuss 
the steady-state regime obtained under constant electri- 
cal excitation. The corresponding set of algebraic equa- 
tions for the steady-state expectation values are reported 
in Appendix [B] In Sec. IIVI the electroluminescence 
spectra are analytically calculated as a function of the 
populations and the appearance of intersubband cavity 
polaritonic resonances in the emission spectra is shown. 
Numerical applications of the theory are presented in 
Sec. |Vl using a specific configuration for the injection 
and extraction electronic reservoirs. The results predict 
the current-voltage characteristics, emission spectra and 
quantum efficiency using different (controllable) param- 
eters for the considered microcavity system. The results 
are critically discussed with respect to the approxima- 
tions of the model. Finally, conclusions and future per- 
spectives are drawn in Sec. IVII 



I. DESCRIPTION OF THE SYSTEM AND 
QUANTUM HAMILTONIAN 

The system under study is described by the following 
second quantization Hamiltonian 



k,(T k,o" k 

k,q,(T k,q,(T 



other- 



(1) 



r 



The energy dispersions of the two quantum well con- 



duction subbands are /kJi(k) 
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-£'12 + §;^i being k the electron in-plane wavevector 
and m* the effective mass (non-parabolicity is here ne- 
glected). The corresponding electron creation fermionic 
operators are c{ ^ j^ and 02^.^., where a is the electron 
spin. ujc{q) = "7= \/Qz + 9^ ^^ ^^^ bare frequency dis- 
persion of a cavity photonic branch as a function of the 
in-plane wavevector q, where c is the light speed, e^ is 
the cavity spacer dielectric constant and Qz is the quan- 
tized photon wavevector along the normal direction, aj, 
is the corresponding photon creation operator, obeying 
bosonic commutation rules. Due to the well-known po- 
larization selection rules of intersubband transitions, we 
omit the photon polarization, which is assumed to be 
Transverse Magnetic (TM). For simplicity, we consider 
only a photonic branch, which is quasi-resonant with the 
intersubband transition, while other cavity photon modes 
are supposed to be off-resonance and can be therefore ne- 
glected in first approximation. The interaction between 



the cavity photon field and the two electronic subbands 
is quantified by the coupling constant 
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where di2 is the intersubband transition dipole along 
the quantum well growth direction, UJ12 = £'12 /fi- is the 
frequency of the intersubband transition, eq the vacuum 
permittivity, L^av is the effective cavity length and S is 
the sample area. For simplicity, we have considered a 
A/2-cavity, with g^ = tt /Lcav being the quantized vec- 
tor along the growth direction. The geometrical factor 

2 
( IL ^ ,V-^+ ^ originates from the TM-polarization nature 

of the transition. Moreover, in Eq. ([2]) we have assumed 
that the active quantum well is located at the antinode 
of the cavity mode field, providing maximum coupling. 
Note that here we have neglected the anti-resonant terms 
of the light-matter interaction and therefore we can de- 
scribe the strong coupling regime for the electrically ex- 
cited system, but not the ultrastrong coupling limit, as 



instead done in Refs. [a, M, W- The Hamiltonian term 
Hother is meant to include all the other interactions: (i) 
electron-phonon interaction; (ii) electron-electron inter- 
action; (iii) electron tunneling coupling to the injection 
and extraction reservoir; (iv) coupling between the cavity 
photon field and the extracavity field. 



II. CLOSED SET OF DYNAMICAL EQUATIONS 
FOR THE ONE-TIME EXPECTATION VALUES 

It is known that due to the cubic light-matter coupling 
term in the Hamiltonian (the product of two fermion op- 
erators and one boson operator) it is not possible to write 
down an exact closed set of equations for the evolution of 
operators, being the Heisenberg equation of motion for 
each product of N operators coupled at least with one 
product of A^ + 1 operators. In other words, the equa- 
tions of motion of the different observables of the system 
form an infinite hierarchy. One approximation method 
that has been used in order to solve this kind of systems 
is the so-called cluster expansion scheme [la. [la. Il7| . It is 
based on a systematic development of expectation values 
of operator products in terms of correlation functions. 

In order to obtain a consistent truncation scheme, 
a pair of fermionic operators has to be considered of 
the same order as a single bosonic operator. In this 
work, we have truncated the hierarchy at the level of 
the product of two excitation operators (i.e., the product 
of four fermion operators). The details of the factor- 
ization are in Appendix [X] The expectation values en- 
tering the present cluster factorization are the electronic 
and photonic populations, the correlation between the 
cavity photon field and the intersubband polarization, 
as well as polarization-polarization correlations. The 
electron occupation numbers in the two quantum well 



conduction subbands 
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are ni^k =< c[^^^^ci^cr,k 

n2,k =< Cj g. i(;C2,(T,k >• Note that, since in the absence of 
a magnetic field all quantities are spin-independent, we 
omit the spin-index in the notation of the averaged quan- 
tities. The cavity photon number is ria^q =< a^Oq >. 
The correlation between the cavity photon field and the 
intersubband electronic polarization is represented by the 
quantity 



Y{q, k) =< alc^ ^C2,a,k+q > 



(3) 



Finally, the polarization-polarization correlation function 
is given by 

X(q+k',k',k) = X! < 4,<T,q+k'Cl,<T,k'cl„,kC2,o-',k+q > • 
a 

(4) 
Note that in the spontaneous photon emission regime, 
y(q, k) can not be factorized: in fact, spontaneous emis- 



sion is incoherent and < Oq >= 0, < c[ ^ k"^2,CT,k+q >= 0, 
meaning that the cavity field and the intersubband polar- 
ization have no definite phase. Loss of coherence due to 
dephasing processes and photonic losses is phenomeno- 
logically quantified by the damping rate Ty- Unlike 
y(q, k), X(k' -I- q, k', k) can be factorized in products of 
non-zero lower-order expectation values of operators. In 
fact, we have X{'k' + q, k', k) = 2n2,k+q(l — "-i,k)^k,k' + 
(5X(k' -I- q, k', k). The first contribution is an uncorre- 
lated plasma term, while 5X(k' -|- q, k', k) describes the 
higher-order correlation, which can be destroyed by de- 
phasing processes quantified by the damping rate Tx- 

The terms in Hother, namely the phonon scattering, 
electron-electron interaction, the coupling to the contact 
reservoirs and the coupling to the external electromag- 
netic field will be treated in an effective way. The car- 
rier non-radiative relaxation (due to phonon-electron and 
electron-electron scattering) is modeled in terms of a sim- 
ple phenomenological relaxation time Tk. Note that the 
role of Coulomb electron-electron interaction on intersub- 
band transitions has been studied, e.g., in Ref. 19. In 
the case of subbands with parallel parabolic dispersion 
(e.g., same effective mass), Coulomb interaction produces 
a moderate renormalization of the intersubband transi- 
tion frequency (^12 and of its oscillator strength, which 
will not be accounted explicitly in the present work. 

Let n1 ^ and n§ j^. be the self-consistent local equilib- 
rium occupation numbers. They are given by Fermi- 
Dirac distributions: 
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where /? = 1/{KT) is the Boltzmann thermal factor, and 
ep is the quantum well self-consistent Fermi level, such 
that: 
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(6) 
The two subbands are coupled to two electronic reser- 
voirs, named respectively left and right contacts. We will 
call F™- [j. the electronic tunneling rate into the k-mode 
of the subband j = 1,2 from the reservoir p = left, right. 
Analogously F™*,^ is defined as the electronic tunneling 
rate from the k-mode of the subband j into the reservoir 
p. The total in-tunncling and out-tunneling rates are 
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The resulting closed system of equations for the one- 
time expectation values reads: 



Left contact 



QW 



Right contact 




In plane wave vector (arb. units) 



FIG. 1: Sketch of the energy dispersion of the two quantum well subbands and of the minibands in the left and right contacts 
in the zero-bias case. Here the system is in thermal equilibrium and the Fermi level in the quantum well is the same as in 
the two contacts. The doping level in the contacts determines the equilibrium density in the quantum well. The subband and 
minibands have an energy dispersion versus the in-plane wavevector, which is a conserved quantity in the planar structure. 
This electronic structure is embedded in a planar microcavity, with a cavity photon mode quasi-resonant to the intersubband 
transition. 
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A. Injection and extraction tunneling rates 



The wave- vector dependent injection and extraction 
rates in Eq. ([7]) can be in principle of different origin. 
Here we give the formal expression for elastic tunnel- 
ing processes conserving the in-plane momentum. Addi- 
tional processes (such as assisted tunneling) can be ac- 



counted for by adding their contribution to the expres- 
sions for Ti\ and T°f to be inserted in Eq. ^. 

As electronic contact reservoirs, we will consider semi- 
conductor doped superlattices, as it is generally the case 
in unipolar quantum cascade devices. 

The chemical potential in each contact is labeled /ip 
with p — left, right. In each reservoir, we will consider 
miniband states with energy E"^^ ^ . In the elastic tun- 



Left contact 
(Extractor) 



QW Right contact 

(Injector) 




In plane wave vector (arb. units) 



FIG. 2: Same as in Fig. [TJ but with an applied voltage bias. Here, the left contact acts as an electronic extractor, while the 
right one is the injector. In the quantum well, non-equilibrium steady-state populations can be established in the two subbands. 



neling process, electron energy and in-plane momentum 
are conserved. The tunneling rate from the contact reser- 
voir into the j-th subband is 
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where Vpj^k^k^ is the tunneling matrix element and k^ is 
in general an index over the electronic states of the mini- 
band with in-plane wave vector k. It can be interpreted 
as the axial electronic wave vector in the case the two 



leads are just bulk contacts. 1/(1 -I- e'^^ p-'^-''^''^"') is the 
Fermi-Dirac occupation number of the electron states in 
the contact. Analogously the tunneling rate from the j-th 
subband of the quantum well into the reservoir p reads 
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is the hole occupation number in the contact. The value 
of r™'"^ can be quantum engineered, depending on the 
specific structure. In particular, by changing the thick- 
ness of the potential barriers, it is possible to tailor con- 
siderably the tunneling matrix element. It is straight- 
forward to see that a simple relationship occurs between 
r™ _k and T;J^^, namely 
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Note that here we have assumed that the bare energy 
dispersion of the electrons in the two subbands is unaf- 
fected. This is valid in the weak light-matter coupling 
regime or when the injector miniband energy width is 



broad enough. For large values of the vacuum Rabi fre- 
quency, the spectral function of the electrons in the sec- 
ond subband is non-trivially modified as well as the tun- 
neling process using a narrow-band injector. This will be 
proved and discussed in detail in a forthcoming paper [lj|. 



III. STEADY-STATE REGIME AND 
OBSERVABLE QUANTITIES 

In this work, we will focus on the steady-state so- 
lutions for the quantities na,q, ^i^k, 'T'2,k, Y{q, k) and 
X(q + q',q',k). Hence, we can set the time derivatives 
equal to zero, transforming the differential system ([7]) 
into an algebraic one. In the steady-state regime, the 
electronic current (number of electrons per unit time) 
through the structure is given by the expression: 

I = ErrkVk-ri'^(i-ni,k) 

k 

= Er2?k(i-"2,k)-r°>2,k. (11) 

k 

The total rate of photons emitted out of the microcavity 
reads 



P = 27^na,q, 



(12) 



where 1/(27) is the escape time of a photon out of the 
microcavity. The quantum efficiency rj is defined as the 
ratio between the photonic current out of the cavity and 
electronic current , i.e., -q = —. 



IV. EMISSION SPECTRA 

In the steady-state regime, the momentum-dependent 
spontaneous photon emission spectra are given by the 



expression: 



In order to determine < a]_(0)aq(i) >, we need to solve 
the following Heisenberg equations of motion 
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q' q' 

Note that here we have omitted the coupling of the electronic injector and extractor reservoirs to the quantity 
< flt (0)cj ^k'^2,cr.k+q >• This coupliug would involve correlations between the quantum well electronic field and the 
contact electronic fields. Since in this paper we are dealing with incoherent electron transport, we will neglect such 
correlations with the contact reservoirs, which are also extremely tricky to tackle. 

Truncating the hierarchy at the level of two excitations (details in Appendix |B]) and taking the unilateral Fourier 
transform {J^ die*'^*) we obtain: 



Sq{t = 0) = na,q = iiu; - tJc(q) + iTsic[))Sq{u;) + 2ix*(q)^q(w) 
Z^{t = 0) = ^r(q,k)=z(w-c^i2+zrz(q))^q(a;)+ix(q)5q(c^)7^. 



(15) 



where S'q(t) =< a^(0)aq(t) >, Zq{t) — J^k < flq(0)c{ ^^ k"^2,o-,k+q > and D represents half the difference between the 
total number of electrons in the fundamental subband and the number in the second one, namely: 
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(16) 



Note that the total density of electrons is 2^k"-i,k + "'2,k, where the 2 factor accounts for the two- fold spin degeneracy 
of the electron states in the conduction subbands. Tg and Tz are phenomenological damping rates for Sq and Zq 
respectively. The analytical solutions are 
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The electroluminescence spectrum is simply 

£q(cj) CX ^Sq{u). (18) 

From the analytical result for Sq{uj), we see immediately 
that emission spectrum is resonant at the two polariton 
frequencies Lo±{q) satisfying the equation 

(w - wi2 + irz){uj - uj,{q) + iTs) - 2x{qfD = 0. (19) 

The quantity ^r = x(q)V2-D is just the vacuum Rabi 
frequency of the present system. At resonance (i.e.. 



^c(q) = ^12), the necessary condition for the appear- 
ance of a strong coupling polaritonic splitting is -D > 



^0 = 



_ (Fs-r. 
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1^, meaning that the total density of elec- 
trons in the fundamental subband must be larger enough 
than the total density in the second. For a vacuum Rabi 
frequency much larger than Tz and Tg, the minimum 
polariton splitting is given by twice the vacuum Rabi 
frequency. 

Note that here the electroluminescence spectral shape 
does not depend explicitly on the spectral properties of 
the injector and extractor reservoirs. The spectrum in 
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FIG. 3: Current density versus applied voltage for different 
values of the intersubband transition energy: E12 ~ 50meV, 
(dashed-dotted line), 100 meV (dashed line) and ISOmeV 
(solid line). Other parameters can be found in the text. 



FIG. 4: Photonic current density versus electronic current 
for different values of the intersubband transition energy: 
E12 = 50meV, (dashed-dotted line), lOOmeV (dashed line) 
and ISOmeV (solid line). Same parameters and range of ap- 
plied voltages as in Fig. |3l 



Eq. P?)l has the same shape as the absorption (in pres- 
ence of the same carrier densities). The dependence on 
the transport is only implicit, being given by the steady- 
state carrier and photon populations. In contrast, in the 
exact solution of the the simplified model of Ref. [Gj, it 
is shown that the electroluminescence spectra is the ab- 
sorption spectrum times the spectral distribution of exci- 
tations in the electronic reservoir, which then acts as an 
electronic filter [l3|. A fermionic approach based on an 
exact diagonalization method 14\ indeed shows that the 
spectral properties of the electronic contact modifies sig- 
nificantly the spectral shape of the electroluminescence 
in the case of narrow band injectors. Hence, the spec- 
trum predicted by Eq. pT|) is valid only for broad band 
injectors. This is not really surprising because, in or- 
der to calculate the tunneling rates, we have considered 
bare electronic states in the quantum well and have only 
considered incoherent population injection and extrac- 
tion processes. 



NUMERICAL APPLICATION 




FIG. 5: Electron occupation number in the fundamental con- 
duction subband as a function of kinetic energy and applied 
voltage. Inset: the integrated density of electrons in the fun- 
damental subband versus voltage. £12 ~ ISOmeV and other 
parameters as in Fig. |3l For eV = E12, the density of elec- 
trons in the first subband is 8.3 x lO^cm"^. 



Here, we apply our theory using realistic parameters 
for a microcavity-embedded quantum cascade electro- 
luminescent source. In order to simplify the algebra, 
we have systematically neglected the photon wavevector 
whenever added to an electronic wavevector. Given the 
huge difference in the typical wavectors of photons and 
electrons, this simplification is safe. Applying this ap- 
proximation, we can obtain a closed set of algebraic equa- 
tion where the variables are the populations in the two 
subbands and in the cavity photonic branch, as shown 
in Appendix [B) This system has been solved numeri- 
cally using a standard Newton method. We achieve nu- 
merical convergence in a relatively fast computation time 



except in the limit of vanishing bias, when the injector 
and extractor are strongly 'misaligned' with the two sub- 
bands. Physically in this case the steady-state situation 
is reached in times very long compared to the dynamics 
of the quantum well system, the photon population is ex- 
tremely small and correspondingly the numerical method 
fails to converge. Anyway this is not a real limitation, 
because we are interested in the behavior of the system in 
presence of a finite voltage bias, producing a significant 
current flow and photonic output. 

In Figs. [Hand [21 we show a sketch of the energy profile 
of the injector and extractor with respect to the quantum 



well subbands respectively without and with an applied 
bias. Specifically, in the numerical calculations we have 
used the following electronic out-tunneling rates: 
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(20) 
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where a = 0.1iJi2, l/F = 0.4ps, £^o,icft and i?o, right are 
the energy offsets of the left and right minibands. The 
in-tunneling rates are determined by applying the rela- 
tion in Eq. (jlOp . In all the simulations, we have taken 

£^0,lcft = £^0, right = 0.5?i(jJi2 and flldt — bright = 3^12- 
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FIG. 6: Same as in Fig. [5] but for the second subband. In- 
set: the integrated density of electrons in the second subband 
versus voltage. For eV = E12, the density of electrons in the 
second subband is 4.3 x lO^cm"^. 
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FIG. 7: Contour plot of the photon occupation (log scale) 
versus the applied voltage and the energy of the bare cavity 
photon mode. 



Note that these are just phenomenological injection 
rates. For the amplitude F, we have considered values 
which are consistent with what realistically obtainable in 
semiconductor intersubband devices. Importantly, in real 
structures F can be considerably quantum engineered 
by changing the barrier thickness and/or the miniband 
structure of the injection superlattices. This is why we 
have not considered a very specific injector configuration 
and taken the simplified expression in Eq. (|20|) with re- 
alistic parameters. 

When a voltage bias is applied, the two reservoirs are 
shifted symmetrically, as shown in Figl^l In all the simu- 
lations, except when otherwise stated we used the realis- 
tic damping parameters Tx = Fy = F5 = Tz = 0.1a;i2, 
7 — 0.05ijJi2, while the temperature is T = 77K. In the 
simulations we have also considered t^ to be independent 
from k and such that - — 0.005wi2, except when oth- 
erwise stated. Note that here we have considered only 
an active quantum well. For quantum cascade struc- 
tures with several active quantum wells repeated in a 
periodic way, the dynamics is similar and the present 
treatment can be generalized without major difhculties. 
In the simulations, the intersubband transition energy 
E12 — huji2 is, except where otherwise stated, equal to 
150meV and the coupling constant x{^) is such that the 
vacuum Rabi frequency is 0.1wi2 for an electron density 
of 5 X lO^^cm"^ (all in the fundamental subband). When 
E12 is changed, the coupling constant is adjusted in or- 
der to keep the ratio between the vacuum Rabi frequency 
and transition frequency constant. The effective mass 
m* has been taken to be one tenth of the bare electronic 
mass. In the numerical calculations, the cavity spacer 
dielectric constant is e^ = 10. For each simulation, the 
resonance in-plane wavector g^es i given by the condition 
^cilres) = '-^12 7 corrcsponds to an internal cavity pho- 
ton propagation angle Ores equal to 70 degrees, where 

tan f^es — Qres / Qz- 

In Fig. [21 we show the current density versus applied 
voltage (between the injector and extractor) for differ- 
ent values of Ei2. The current-voltage profile is char- 
acteristic of an unipolar quantum cascade light emitting 
diode. The current grows superlinearly in the voltage re- 
gion where the injector Fermi level approaches the second 
subband. The current is bigger for smaller E12 because, 
keeping the injection rate F constant (but all the inter- 
nal rates of the system proportional to E12), the injection 
and extraction processes become the dominant processes. 
Note that an increase of the nonradiative relaxation rate 
I/t produces a nearly proportional increase of the elec- 
tronic current (not shown) . The rates of emitted photons 
per unit area (integrated all over the in-plane wavevec- 
tors) are shown in Fig. 2] as a function of the flowing 
current, showing an approximately linear behavior. 

Fig. Elandinishow contour plots of the electron occupa- 
tion numbers of the first and second subband respectively 
as a function of the applied voltage and of the kinetic 
energy. The insets in Fig. [5] and [5] show respectively 
the integrated density of electrons in the first and sec- 



ond subband. It is apparent that with increasing voltage 
the population in the first subband decreases, while the 
population in the second subband increases. 

When the injector Fermi level becomes aligned with 
the second subband, as expected, the carrier occupation 
numbers in the two subbands are considerably out of 
equilibrium. The decrease of the first subband carrier 
occupation numbers is beneficial for the radiative effi- 
ciency of the spontaneous emission, because the influence 
of Pauli blocking is reduced. Moreover, in the considered 
conditions, the density of electrons in the ffi'st subband 
is still considerably larger than in the second subband, 
thus producing a large vacuum Rabi coupling and effi- 
cient emission rate. 

Fig. [7] contains a contour plot of the cavity photon 
occupation number versus the bare photon energy, show- 
ing that the maximum of emission is obtained when the 
bare photon energy is resonant with the intersubband 
transition, as expected and as observed experimentally 
[l3,[l3l- With the considered parameters, the density of 
electrons in the first subband is high enough to be in the 
strong coupling regime, as depicted in Fig. [8l where the 
anticrossing of two polariton branches is clearly present. 
The minimum polariton splitting, given by the expression 
2xiq)V2D is reported in Fig. [9] as a function of the ap- 
plied bias. With increasing voltage, the population differ- 
ence D = X^k ^k = X^k'^i.k ~ '^2,k diminishes. This re- 
sults in a decrease of the vacuum Rabi frequency and con- 
sequently of the polariton splitting. This high-excitation 
feature has been already observed in experiments [13|,ll8| 
and can not be described within a bosonic approach[y|, 
which can be applied only to the low excitation density 
case [205. 

It is interesting to analyze the quantum efficiency 77, 
defined as the ratio between the photonic emission rate 
and the electronic current, namely 



»7 = 
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Ekrrk^i^k - r- (1 - ni,k) ■ 



(21) 



In Fig. 1101 we plot the quantum efficiency rj at 
eV — E12 versus the vacuum Rabi frequency ^Ir at the 
same voltage (log- log scale). In the simulations, the vac- 
uum Rabi frequency has been varied by changing the cou- 
pling constant x(q)- In ^i realistic quantum engineered 
device, x(q) can be tailored in different ways. For ex- 
ample, by growing the active quantum wells in a spatial 
region where the cavity mode field is very small, it is 
possible to quench dramatically the value of x(q). More- 
over, by using different shape of quantum wells, it is also 
possible to tailor the transition dipole di2- Fig. [TOl shows 
that in the weak coupling regime (small values of Qri.) the 
efficiency grows like fl^. In the strong coupling regime, 
the efficiency becomes impressive and then tends to satu- 
rate. It is apparent that the radiative efficiency smoothly 
increases passing from the weak to the strong coupling 
regime. This crossover occurs because the radiative effi- 
ciency depends on the spectrally integrated emission and 




CO (q) (in units of co^ 

c ^ 1 



FIG. 8: Contour plot of the electroluminescence (arb. units) 
as a function of the bare cavity photon energy uJc{q) and of 
the emission frequency u for an applied voltage eV = 0.5-Bi2. 
The anticrossing of the two intersubband polariton branches 
is apparent in electroluminescence spectra. 



it is therefore insensitive to the sudden appearence of the 
polariton doublet in the strong coupling emission spectra. 

This results are in qualitative agreement with the ana- 
lytical solutions of the simplified model in Ref. 6, where 
only the bright intersubband states are considered and 
where the electronic reservoir is modeled with a bath of 
harmonic oscillators. As shown in Fig. [TTl the nonradia- 
tive population relaxation rate 1/r has the most signifi- 
cant effect. In the considered regime of parameters, the 
efficiency is proportional to r. 

It is interesting to compare our results for this mi- 
crocavity system with the standard free space case. In 
the free-space case, the photon current, obtained by 
applying the Fermi golden rule, is given by the for- 
mula P = ^%^Ti^Ek"-2,k(l - »^i,k)- As it is well 
known, the free-space radiative efficiency dramatically 
decreases with the intersubband emission wavelength due 
to the (^12^12 dependence of the spontaneous emission 
rate (d^2 oc 1/cl;i2, so the spontaneous emission rate 
scales effectively as (^12). In the mid-infrared, by using 
the same parameters, for a transition of 150 meV, the 
quantum efficiency is of the order of 10""' — 10~^. Hence, 
it is clear from our results that a strong coupling light- 
emitting diode based on a planar microcavity system can 
provide a dramatic enhancement with respect to the free 
space case (even three orders of magnitude for the larger 
vacuum Rabi frequency case). 



VI. CONCLUSIONS AND PERSPECTIVES 

In conclusion, we have presented a quantum theoreti- 
cal study of the quantum well intersubband electrohimi- 
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FIG. 9: Minimum polariton splitting as a function of the 
applied voltage. E12 = 150meV and other parameters can be 
found in the text. 
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FIG. 10: Quantum efficiency versus the corresponding vac- 
uum Rabi frequency at the voltage eV = E12 . The three lines 
are obtained with different values of the coherence damping 
coefficients: Fx ~ Fy = 0.1lji2 (solid line), 0.05(x;i2 (dashed 
line) and 0.025a;i2 (dashed-dotted line). E12 = ISOmeV. 



nescence from a semiconductor microcavity in the inco- 
herent transport regime, i.e. when the couphng to the 
electronic contacts concerns only the electron popula- 
tions. The problem has been tackled starting from the 
fermionic electron Hamiltonian for the two-subbands and 
by using a cluster factorization method to truncate the 
infinite hierarchy of dynamical equations for the relevant 
expectation values of operator products. At the present 
level of approximation, we have been able to describe 
the incoherent electron transport through the dynamics 




10 10' 

CI (in units of Q) ) 



10 



FIG. 11: Quantum efficiency versus the corresponding vac- 
uum Rabi frequency at the voltage eV — E\2- The three 
lines are obtained with different values of the non-radiative 
relaxation rate: - — 0.01lji2 (dashed line), 0.005lji2 (solid 
line) and 0.0025a;i2 (dashed-dotted fine). E12 — 150meV. 



of the electronic subband populations and the electrolu- 
minescence through the dynamics of the cavity photon 
population, which is coupled to the correlations between 
the electromagnetic field and the intcrsubband polariza- 
tion. We have discussed the limits of applicability of 
the present approach, which neglects the impact of the 
vacuum Rabi coupling on the quantum well carrier spec- 
tral function and any correlation between the quantum 
well and the contact reservoirs. The analogies and differ- 
ences with the exact predictions of the simplified model 
in Ref. lo have been critically and extensively discussed. 
We have shown the appearance of cavity polariton reso- 
nances in the emission spectra, when a large density of 
electrons occupies the fundamental subband. We have 
described how the vacuum Rabi splitting decreases with 
increasing voltage and described the photonic output in 
the different transport conditions. Our results show that 
even in presence of non-radiative relaxation and Pauli 
blocking, the quantum efficiency of the microcavity in- 
tcrsubband electroluminescence can be considerably en- 
hanced by increasing the vacuum Rabi frequency. A more 
refined treatment |14| based on a fermionic exact diago- 
nalization method shows that under certain conditions 
the strong vacuum Rabi coupling regime affect consider- 
ably not only the dynamics of the intcrsubband polariza- 
tion (hence the absorption spectrum), but also the quan- 
tum well electron spectral properties and consequently 
the tunneling transport using narrow-band injectors. As 
future perspective, this could be exploited to further im- 
prove the quantum efhciency of microcavity intcrsubband 
emitters and for the eventual realization of intcrsubband 
polariton lasers. 
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As stated in the main body of the paper we used a cluster expansion and truncation scheme to obtain a closed set 
of equations. Here we briefly review the principles of this method following [1^, \1Q, [l3| and apply it to the actual 
case. 

If we consider each bosonic operator or each pair of fermionic operators as an excitation operator and we write the 
expectation value of an N excitation operator as < N >, then the Heisenberg equation of motion takes the form: 

i—<N>^ T[< iV >] + y[< iV + 1 >] 

where the N-excitation expectation value is coupled to higher order quantities via the functional V. An N-excitation 
truncation scheme is obtained if we factorize all the expectation values of more than N excitation in all the possible 
ways and considering the sign exchange for the fermionic operators in order to obtain a factorized quantity that 
respects the commutation and anticommutation properties of the original quantity. 

We are interested in incoherent emission only, so the only nonzero one excitation operators we consider are < 
c[ ^ k'-i.o-.k > ^-nd < C2 g. k'^2,cr,k >• We factorized the 3 excitations operators in the following way: 

< aqCl,<T,k4,<T,k'4,a',k"C2,^',k"' > = - < aq4,^,k'Cl,<T,k >< 4,o-',k"C2,o-',k"' > + < aq4,o-',k"Cl,^,k >< 4,a,k'C2,o-',k"' > 

= - < aq4,o^,k'Cl,(T,k > '5k",k""^2,k"+ < «q4,<T,k"Cl,(T,k > '5k%k"'5<T,o-'rt2,k' 

< '^q4,cr,kCl,cr,k'Cl,o-',k"4,o-',k"' > = ^ < "q4,o-,kCl,(T',k" >< Cl,o-,k'4,o-',k"' > + < «q4,o-,kCl,(T,k' >< Ci,o-',k"4,o-',k"' > 

= - < aq4,o^,kCi,<T,k" > '5k',k""5cr,o^'(l - ni,k')+ < aq4,<T,kCi,<T,k' > Sk"M"'{^ - "i.k") 
For the two-time quantities in the calculation of luminescence spectrum, we proceed analogously and obtain: 

< aq(0)aq'4,a,k+qC2,<T,k+q' > = < a]j(0)aq' >< 4,^,k+qC2.<T,k+q' > ^q,q' , 

< aq(0)aq/c}^^^Cl^<j,k+q-q' > = < a]j(0)aq/ >< 4_^,j.Cl,cr,k+q-q' > (^q^q' • 

APPENDIX B: ALGEBRAIC EQUATIONS FOR THE STEADY-STATE REGIME 

In the steady-state regime, neglecting the photonic wavevector into sums over electronic wavevectors, the system 
of equations ^ reduces to the following system of algebraic equations. 

= I Sq(7 + Tx) + (Jy + 2Dx(q)^ '^ j ""''' "D ^ " '' n. ^i^k^i^k - T^^{1 - ni,k)) ■. 

'^ = ^l^~r^T — i^-Dk) + o)( + ri,k'^i,k-rik(l-'^i,k)) 

q Gqlx 2 T-k 

ixly ^ trq V 'I 

= ■ —-ntn2M + n\{^-n2M)-— ^-rr^ni,k-fri\(l-ni,k), 

Tk Tk ' ' 



where 
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-Fie = n2,kil - ni,k) 
D = J2Dk 



F 



'l.k 



"■2,k 



Gc 



E^^ 



Ei? is calculated by inverting the relation 
y^»-i,k + n2,k 



27r?i2 



1 



Bq — Ty 



exp/?(u;i(k) ~ ep) + 1 

1 
exp/3(a;2(k) - e^) + 1 
2x(q)' 



-D 



Wc(q) - Wi2 

f f \ \2 , /-r , 2x(q)^^ ^2 
(wc(q) - W12) + (ry H ) . 



de 



exp/3(e - ef) + 1 exp/3(e + -^12 - ef) + 1' 



Discretizing the electronic and photonic wavevectors on a grid of respectively iVk and Nq points, we have a system 
of 2A'k + Nq equations that can be numerically solved, e.g., with a Newton algorithm. 



[1] Intersubband Transitions in Quantum Wells: Physics 

and Device Applications I, edited by H. C. Liu and F. 

Capasso, Semiconductors and Semimetals Vol. 62 (Aca- [11] 

demic Press, San Diego, 2000). 
[2] D. Dini, R. Kohler, A. Tredicucci, G. Biasiol, and L. [12] 

Sorba, Phys. Rev. Lett. 90, 116401 (2003). 
[3] A. A. Anappara, A. Tredicucci, G. Biasiol, L. Sorba, 

Appl. Phys. Lett. 87, 051105 (2005). [13] 

[4] A. A. Anappara, A. Tredicucci, F. Beltram, G. Biasiol, 

L. Sorba, Appl. Phys. Lett. 89, 171109 (2006). 
[5] C. Ciuti, G. Bastard, L Carusotto, Phys. Rev. B 72, [14] 

115303 (2005). [15] 

[6] C. Ciuti, L Carusotto, Phys. Rev. A 74, 033811 (2006). [16] 

[7] S. De Liberato, C. Ciuti, L Carusotto, Phys. Rev. Lett. 

98, 103602 (2007). [17] 

[8] J. Faist, F. Capasso, D.L. Sivco, C. Sirtori, A. L. 

Hutchinson, A.Y. Cho, Science 264, 5553 (1994). [18] 

[9] R. Kohler, A. Tredicucci, F. Behram, H.E. Beere, E.H. [19] 

Linfield, A.G. Davies, D.A. Ritchie, R.C. lotti, F. Rossi, 

Nature 417, 156 (2002). [20] 

[10] R. Colombelli, K. Srinivasan, M. Troccoli, O. Painter, C. 



F. GmachI, D. M. Tennant, A. M. Sergent, D. L. Sivco, 
A. Y. Cho, F. Capasso, Science 302, 1374 (2003). 
R. Colombelli, C. Ciuti, Y. Chassagneux, C. Sirtori, 
Semicond. Sci. Technol. 20, 985 (2005). 
L. Sapienza, A. Vasanelli, C. Ciuti, C. Manquest, C. Sir- 
tori, R. Colombelli, and U. Gennser, Appl. Phys. Lett. 
90, 201101 (2007). 

L. Sapienza, A. Vasanelli, R. Colombelli, C. Ciuti, 
Y. Chassagneux, C. Manquest, U. Gennser, C. Sirtori, 

Phys. Rev. Lett. 100, 136806 (2008)^ 

S. De Liberato and C. Ciuti, arXiv: 08023091] 

J. Fricke, Ann. Phys. 252, 479 (1996). 

J. Fricke, V. Meden, C.Wohler, K. Schonhammer, Ann. 

Phys. 253, 177 (1997). 

M. Kira, W. Hoyer, T. Stroucken, S. W. Koch, Phys. 

Rev. Lett. 87, 176401 (2001). 

L. Sapienza, (Ph. D. thesis. University of Paris 7, 2007). 

D. E. Nikonov, A. Imamoglu, L. V. Butov, and H. 

Schmidt, Phys. Rev. Lett. 79, 4633 (1997). 

M. F. Pereira, Jr. , Phys. Rev. B 75, 195301, (2007). 



